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Abstract 

We perform a coupled-channel analysis of pseudodata for the D'^ —)• Dalitz plot. The 

pseudodata are generated from the isobar model of the E791 Collaboration, and are reasonably re¬ 
alistic. We demonstrate that it is feasible to analyze the high-quality data within a coupled-channel 
framework that describes the final state interaction of D~^ —)• as multiple rescatterings 

of three pseudoscalar mesons through two-pseudoscalar-meson interactions in accordance with the 
two-body and three-body unitarity. The two-pseudoscalar-meson interactions are designed to re¬ 
produce empirical tt-tt and ttK scattering amplitudes. Furthermore, we also include mechanisms 
that are beyond simple iterations of the two-body interactions, i.e., a three-meson force, derived 
from the hidden local symmetry model. A picture of hadronic dynamics in ^ de¬ 

scribed by our coupled-channel model is found to be quite different from those of the previous 
isobar-type analyses. For example, we find that the —)• A'^vr+Tr"’' decay width can get almost 
triplicated when the rescattering mechanisms are turned on. Among the rescattering mechanisms, 
those associated with the p{n0i)K^ channel, which contribute to D'^ —)■ only through a 

channel coupling, give a large contribution, and significantly improve the quality of the fits. The 
s-wave amplitude from our analysis is reasonably consistent with those extracted from the 
E791 model independent partial-wave analysis; the hadronic rescattering and the coupling to the 
/9(770)A*^ channel play a major role here. We also find that the dressed decay vertices have 
phases, induced by the strong rescatterings, that strongly depend on the momenta of the final 
pseudoscalar mesons. Although the conventional isobar-type analyses have assumed the phases to 
be constant, this common assumption is not supported by our more microscopic viewpoint. 
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I. INTRODUCTION 


With the advent of charm and B-factories, a large amount of data for charmed-meson 
decays have been accumulated in the last decades. Among a number of physical interests, one 
appealing aspect of studying these charmed-meson decays is that we can gain information 
about interactions between light mesons and resonances. This was particularly highlighted 
by the E791 Collaboration’s report on their identihcation of the a meson in the Dalitz 
plot of the —)■ 7r“7r’''7r’'' decay [1]. A similar analysis was also made for the —)■ 
A“7r’''7r+ decay to identify the k resonance [2]. These Endings triggered further analyses 
of the —)■ K~7r~^7r~^ Dalitz plot data, paying special attention to the K~tt^ s-wave 
amplitude, as follows. Oiler [3] analyzed the E791 data [2] using the Kti 1=1/2 (/: total 
isospin) s-wave amplitude based on the chiral unitary approach [4], instead of Breit-Wigner 
functions for the K~tt^ s-wave used in the E791 analysis, and obtained a reasonable fit. 
The 1=3/2 s-wave was not considered in his analysis. The E791 Collaboration reanalyzed 
their data without assuming a particular functional form for the K~7r^ s-wave amplitude. 
Rather, they determined it bin by bin, which they call model independent partial-wave 
analysis (MIPWA) [5]. An interesting finding in the MIPWA was that the obtained JP“7r+ 
s-wave amplitude has the phase that depends on the K~ti^ energy in a manner significantly 
different from what is expected from the Watson theorem combined with the LASS empirical 
amplitude [6], assuming the J=l/2 dominance. Edera et ah suggested that this difference 
can be understood as a substantial mixture of the Kn 1=1/2 and 3/2 s-wave amplitudes [7]. 
This idea was implemented in the Dalitz plot analysis done by the FOCUS Collaboration [8]. 
They parametrized the K~7r^ s-wave amplitude in terms of the iC-matrix of 1=1/2 and 
3/2 that had been fitted to the LASS amplitudes [6, 9]. They found that the D+ —)■ 
A“7r+7r+ Dalitz plot can be well fitted with the K~7r~^ s-wave amplitude in which the 
1=1/2 and 3/2 components interfere with each other in a rather destructive manner. The 
FOCUS Collaboration also has done a MIPWA in a subsequent work [10] to find results 
similar to those of the E791 MIPWA. The quasi-MIPWA has also been done by the CLEO 
Collaboration [11]. Their new finding was that 1=2 tttt nonresonant amplitude can give 
a non-negligible contribution. Meanwhile, an analogous decay, D+ —)■ has been 

analyzed by the BESIII Collaboration [12]. An interesting finding was that the p(770)A^ 
channel gives by far the dominant contribution. This implies that, although the p(770)A^ 
contribution is not directly observed in the D~^ —)■ K~7r~^7r~^ Dalitz plot, it can play a 
substantial role in the hadronic final state interactions (ESI) of the —)■ K~7r~^7r~^ decay, 

considering that the two decays share the same hadronic dynamics to a large extent. 

Many of the previous analyses of the D~^ —)■ K~7r^'ir^ Dalitz plot have been done with 
the so-called isobar model in which a D”'"-meson decays into an excited state R {k, K*, K^-, 
etc.) and a pseudoscalar meson. The R subsequently decays into a pair of light pseudoscalar 
mesons, while the third pseudoscalar meson is treated as a spectator. The propagation of R 
is commonly parametrized by a Breit-Wigner function. The total decay amplitude is given 
by a coherent sum of these isobar amplitudes supplemented by a fiat interfering background. 
On the other hand, as mentioned in the above paragraph, some analysis groups modified 
the conventional isobar model to use the iP-matrix or chiral unitary model for the iP“7r+ 
s-wave amplitude so that the consistency with the LASS data for K~tt^ —)■ and 

with the two-body unitarity is maintained. Meanwhile, in (Q)MIPWA, the K~7r^ s-wave 
amplitude is solely determined by the D~^ —)■ Dalitz plot data. We will refer 

to these models, which do not explicitly consider three-meson-rescattering required by the 


2 


three-body unitarity, as isobar-type models. The basic assumptions common to all of these 
models are that the spectator pseudoscalar meson interacts with the others very weakly, 
and/or —)■ Rc (c: spectator meson) vertices with complex coupling constants can absorb 
such effects. 

Although each analysis group has obtained a reasonable £t to its own data, their results 
are not necessarily in conformity with theoretical expectations. For example, the E791 
Collaboration reported that the phase of the 1=1/2 Ktt p-wave amplitude used in their 
MIPWA is, according to the Watson theorem, not consistent with that of the LASS analysis 
in the elastic region [5]. This may be originated from either or both of two possible reasons 
below, each of which signals a serious problem in the basic assumptions underlying in the 
analysis model. One possible reason is that the 1=1/2 Kir p-wave amplitude in the E791 
analysis is given by a coherent sum of Breit-Wigner functions for .A* (892) and A* (1680), 
and it is not consistent with the LASS data to the required precision. Or the amplitude 
does not satisfy the two-body unitarity not only formally but also quantitatively. Another 
possibility is that the (neglected) rescattering of the spectator meson with the other mesons 
actually plays a substantial role to generate an energy-dependent phase so that the Watson 
theorem does not hold in reality, and the analysis model tried to fit it. Not only the E791 
analysis but also all the previous analyses mentioned above would share the same problem, 
and this seems to indicate a need for going beyond the conventional isobar-type analysis; a 
unitary coupled-channel approach. In order to extract from data a right physics, e.g., K~ir^ 
s-wave amplitude from K~ir^ir^ Dalitz plot data, one needs to use a theoretically 

sound model so that unnecessary model artifacts do not come into play. 

Recently, we have developed a unitary coupled-channel framework for describing a heavy- 
meson decay into three light mesons [13]; both the two-body and three-body unitarity are 
maintained. In the reference, we studied the extent to which the isobar-type description of 
heavy-meson (or excited meson) decays is valid by analyzing simple pseudodata. We found 
a signihcant effect of the channel couplings and multiple rescattering on the Dalitz plot 
distributions. This study has been extended to an analysis of pseudodata for excited meson 
photoproductions [14]. In the present work, we apply the formalism of Ref. [13] with some 
modifications to a realistic case. Thus, we will perform a coupled-channel analysis of the 
—)■ A“7r+7r+ Dalitz plot pseudodata generated from the E791 isobar model [5]. To the 
best of our knowledge, this is the first coupled-channel Dalitz plot analysis of a D-meson 
decay into three pseudoscalar mesons. We will demonstrate that a quantitative coupled- 
channel partial-wave analysis of the —)■ A“7r’''7r’'' Dalitz plot is feasible. Then we will 

examine the hadronic dynamics in the ESI of the —)■ A“7r+7r+ decay within the coupled- 

channel model. We will study how the partial-wave amplitudes and their fit fractions are 
different between an isobar-type model and a model that includes the three-body scattering. 
We also examine contributions from the rescattering and channel couplings to the Dalitz plot 
distribution. Through these investigations, we address the validity of the above-mentioned 
basic assumptions of the isobar-type model from this more microscopic viewpoint. 

So far, the three-body FSI for the —)■ A“7r+7r+ decay has been explored by Magalhaes 

et ah [15] (see also Guimaraes et ah [16]). They were concerned with whether the difference 
in the phase between the A“7r+ s-wave amplitudes from MIPWA (E791 [5] and FOCUS [10]) 
and the Kir 1=1/2 s-wave amplitude from the LASS analysis can be understood as a result of 
the FSI. They calculated the K~ir^ s-wave amplitude in the D'^ —)■ K~ir^ir^ decay with only 
the Kir 1=1/2 s-wave scattering amplitude of the chiral unitary model that had been fitted to 
the LASS amplitude in the elastic region; other partial waves as well as inelasticities were not 
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taken into account. An interesting finding in their work was that the K~7t~^ s-wave amplitude 
originated from the weak vector current and the subsequent s-wave rescattering is in a fairly 
good agreement with the phases from the MIPWA in the elastic region; the qualitative 
feature of the modulus of the MIPWA amplitude is also described. This hnding was further 
conhrmed in their subsequent works [17] where p{770) contribution was also partly taken 
into account. Although their results are interesting and suggestive, they should be looked 
with a caution for the reasons below. First, their calculated amplitudes are still qualitatively 
in agreement with those from the MIPWA in the elastic region, and more rehnements are 
needed particularly in the inelastic region. There will be a delicate interplay between new 
mechanisms to be included for the rehnements and the already existing mechanisms, and 
thus it is not clear if their current hndings persist after the rehnements. Second, the authors 
of Refs. [15, 17] treated the K~tt^ s-wave amplitudes from the MIPWA as data, and did not 
analyze the Dalitz plot directly. However, we note that the MIPWA s-wave amplitude was 
obtained under the assumption that all the other partial-wave amplitudes are basically not 
changed by the FSI. In case the s-wave amplitude is signihcantly modihed by the FSI, the 
other partial-wave amplitudes are also probably modihed, and the MIPWA s-wave amplitude 
may no longer be compatible with them. In order to fully examine the FSI ehects on each 
of the partial-wave amplitudes for the —)■ decay, it is hrst necessary to analyze 

the Dalitz plot with a model that takes account of all relevant partial waves and the FSI, 
thereby extracting the partial-wave amplitudes, as we will do in this work. Then, the FSI 
ehects on the extracted amplitudes can be studied. 

In our analysis, we use two-pseudoscalar-meson interactions that generate unitary ampli¬ 
tudes for Kti and tttt scatterings, and consider resonances (k, K*, p, etc.) as poles in the am¬ 
plitudes. The two-body interactions are htted to the LASS and CERN-Munich [18-20] data. 
With the two-body interactions, we solve the Faddeev equation for a three-pseudoscalar- 
meson scattering to obtain an amplitude that respects the three-body unitarity and thus 
channel couplings. This amplitude is used to describe the FSI of the decay. The pseu¬ 
dodata are htted by adjusting the strengths and phases of (bare) —)■ Rc vertices, while 

the two-pseudoscalar-meson interactions are hxed as those obtained by htting the two-body 
scattering data. In this way, we will examine the extent to which we can ht the D^-decay 
pseudodata, keeping the consistency with the two-body scattering data for all partial waves 
considered in our model. This is in contrast with the previous D~^ —)■ K~7r~^7r~^ analyses 
where some of the resonance parameters were also adjusted along with the —)■ Rc ver¬ 

tices. We consider the 1=1/2 Ktt s-, p-, and d-waves as commonly included in the previous 
analyses. We also explicitly include the 1=3/2 Kir s-wave {1=2 tttt s-wave) as has been 
done so in the FOCUS [8] (CLEO [11]) analysis. We also consider the 1=1 tttt p-wave where 
the p(770) plays a major role. This partial wave has not been considered in the previous 
Dalitz plot analyses because it does not directly decay into the K~TT~^Tr~^ hnal state. How¬ 
ever, this partial wave can still contribute to —)■ K~7r~^7r~^ through the channel couplings. 

Considering the BESHI analysis mentioned above, the p(770)JF channel is expected to play 
a substantial role also here, and we will see that this is indeed the case, at least within our 
analysis. 

Now we discuss the last piece of our model. The FSI of the —)■ K~TT~^Tr~^ decay 

is a three-pseudoscalar-meson scattering. Generally in a three-body scattering, there can 
exist a mechanism that cannot be expressed by a combination of two-body mechanisms, 
i.e., a three-body force. Based on the hidden local symmetry (HLS) model [21], in which 
vector and pseudoscalar mesons are implemented together in a chiral Lagrangian, we can 
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actually derive a “three-meson force” essentially in a parameter-free fashion, up to form 
factors we include. Thus we consider some of the HLS-based three-meson force acting on 
important channels, and examine how they play a role in the —)■ decay. If two- 

pseudoscalar-meson interactions are well determined by precise two-body scattering data, 
the D —)■ Ktttt decays and also other decay modes such as D ^ tititi could serve as a ground 
to study the three-meson force. 

The organization of this paper is as follows. In Sec. II, we discuss our coupled-channel 
model and present formulas to calculate the —)■ decay amplitude. Then in 

Sec. Ill, we present numerical results from our analyses of the two-pseudoscalar-meson scat¬ 
tering data and of the —)■ it'“7r+7r+ decay Dalitz plot pseudodata. Finally, we give a 
summary and future prospects in Sec. IV. A derivation and the resulting expressions for the 
three-meson force, and model parameters are presented in the Appendixes. 


II. FORMULATION 

We have already developed a formalism to describe a heavy-meson decay into three light 
mesons in Ref. [13]. In the present work, we basically use the same formalism with some 
modihcations. Thus, here we just present expressions that are needed in the following 
discussions, specifying the modifications we make for this work. For derivations of most 
of the expressions, see Sec. II of Ref. [13]. Our formalism can be regarded as a three- 
dimensional reduction of a fully relativistic formulation. Because we deal with scatterings 
of light particles, i.e. pions, one may question the validity of this approximation. Although 
a legitimate concern, we made an argument on this along with improvements needed in 
the future in Sec. V of Ref. [13]. In our formalism, we first construct a two-light-meson 
(tt, K) interaction model that is subsequently applied to three-light-meson scattering. The 
following presentations are also given in this order. 


A. Two-light-meson scattering model 

We describe two-light-meson scatterings with a unitary coupled-channel model. For ex¬ 
ample, we consider tttt and KK channels for a tttt scattering, and tiK and rj'K channels 
for 1=1/2 s-wave tiK scattering. We will specify channels for each partial wave later. We 
model the two-light-meson interactions with resonance(i?)-excitation mechanisms or contact 
interactions or both. We choose this parametrization of the two-meson interactions so that 
it can be easily applied to a three-meson scattering. A difference from Ref. [13] is that we 
include the contact interactions here but we did not in Ref. [13]. This is because we include 
partial-wave amplitudes that do not have resonances in analyzing the —)■ decay 

Dalitz plot. Besides, we use a parametrization for R ^ ab (a, b: pseudoscalar mesons) vertex 
functions that are different from those used in the previous work and, with this parametriza¬ 
tion, we need to include contact interactions in addition to i?-excitation mechanisms in order 
to obtain reasonable hts to some empirical partial-wave amplitudes. 

We consider a partial-wave ab —?■ a'b' scattering (see Fig. 1 for a diagrammatic repre¬ 
sentation) with total energy E, total angular momentum L, total isospin I. We will also 
denote the incoming and outgoing momenta by q an q', respectively; q = \q\ throughout this 
paper, except for the Appendix. When q is on-shell, it is related to E hy E = Ea{q) + Eh{q) 
and Ea{q) = ma being the mass of a. First let us consider a ab ^ a'b' scat- 
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FIG. 1. Diagrammatic representation of two-light-meson partial-wave scattering, ab —^ a'b'. (a) 
two-meson interaction potentials; (b) two-meson partial-wave scattering amplitude; (c) two-meson 
partial-wave scattering amplitude from separable contact interactions; (d) dressed R ^ ab decay 
vertex; (e) the inverse of dressed R Green function. 


tering that is not accompanied by a resonance excitation. In this case, we use a separable 
two-light-meson interaction potential [hrst term of rhs of Fig. 1(a)] as follows: 


Va'b',abW,(l) = W%{(i)K'V^^b WabiQ), 


( 1 ) 


where is a coupling constant, and w^^{q) is a vertex function. We use the following 

parametrization for w^^{q): 


Wab (q) 




1 2 + L /2 


1 + 



( 2 ) 


where is a cutoff parameter. Meanwhile, we allow an exception for L=0, 1=2 tttt scat¬ 
tering for which we use a different parametrization for the vertex function: 




1 -h h'{q/mnY 
En{q) 


1 

.1 + 


2 


3 


( 3 ) 


where an additional coupling constant h' has been introduced to obtain a reasonable £t to 
data. For a later convenience, let us define by 


Wal (?) 


—=w^b (?) (if ® ^^cl b are identical particles), 

v2 

Eabi^) (otherwise). 


( 4 ) 


With the above interaction potential, the partial-wave amplitude is given as follows [see 
Fig. 1(c) for a diagrammatic representation]: 

ta'b',abiQ'^ rE) = Y^ (q), (5) 

a"b" 
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LI 

^a'b',ab 



dq 



E - Ea{q) - Ek{q) + ie ' 


( 6 ) 

( 7 ) 


The amplitude of Eq. (5) can contain resonance pole(s), in general. 

Now we extend the model to include explicit degrees of freedom for excitations of spin-L 
isospin-J resonances that contribute to the ab —)■ a'b' scattering. In this case, a two-light- 
meson interaction potential includes bare i?-excitation mechanisms in addition to the contact 
potential defined in Eq. (1) as follows [see also Fig. 1(a)]: 


Vlb',abW, E) = q) + Y^ fIIb',RW) 


1 


R 


E -rriR 


fR,abi(l) 


( 8 ) 


where ttir is the bare mass of R] fabjiiq) denotes a bare R ^ ab vertex function and 

fklabi^) = fkb%<l)- Following Ref. [13], we introduce that is related to f^h^iq) as 

follows: 



J_JLI 

flbAQ) 


ab,R 


(g) 


(if a and b are identical particles), 
(otherwise). 


( 9 ) 


Then we employ the following parametrization for the bare vertex function: 


fLI 
J ab,R 


(g) = gab,R 


m-n 


^/rriLEJkq)E^ [l -h (q/cab,R 


l+(L/2) 



( 10 ) 


where gab,R and Cab,R are coupling constant and cutoff, respectively. This parametrization is 
different from that used in Eq. (35) of Ref. [13], and has a proper kinematical factor. With 
the potential in Eq. (8), the resulting scattering amplitude is given by [Fig. 1(b)]: 


Tkb'^ab{g 1 Q'l E) — t^k',ab{g 1 Q'i E) -I- ^ fkb',R'{Q'i E)Tk/^R{E)fk^^^{q', E) 

R',R 


( 11 ) 


where the first term is the scattering^amplitude from the contact interactions only, as has 
been defined in Eq. (5). The symbol fab,R denotes the dressed vertex that describes the bare 
R ^ ab decay followed by ab rescattering through the contact interactions. Expressions for 
fab,R and fn^ab are [Fig. 1(d)]: 


JlbAg] E) = f^bUiq) + J2 


a'b' 


tala'b'ig,g'',E) f^i',Rig') 

E - Ea'{q') - Eb'{q') + ie ’ 


fkUig-, E) = fk!abig) + 

a'b' 



fkk'b'W) ta'i',abiQ'^g'^E) 

E - Ea'{q') - Eb>{q') + ie ' 


( 12 ) 

(13) 


In Eq. (11), the dressed Green function for i?, ^{E), has been introduced, and is given 

by [Fig. 1(e)]: 

[(r^^(E))-i]^,^ = {E- mn)6n',R - E^f^(E), (14) 
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where is the self-energy of i?, and is defined by 


ab 



q^dq 


E - Ea{q) - Eb{q) -f ie' 


(15) 


In case ab a'b' interaction is given by only resonant mechanisms [no first term in Eq. (8)], 
which is the case in Ref. [13], the corresponding scattering amplitude is obtained from 
Eq. (11) by dropping the first term, and replacing the dressed vertex (/) with the bare one 
(/) in Eqs. (11) and (15). 

The partial-wave amplitude, Ty^, in Eq. (11), is related to the R-matrix by 

SaiabiE) = Vli = 1 - 2mpab Tyabido, do] E) , (16) 

where qo is the on-shell momentum that satisfies E = Ea{qo)+Eb{qo), and pab = qoEa{qo)Eb{qo)/E 
is the phase-space factor. The phase shift and inelasticity are denoted by 5li and pu, re¬ 
spectively. 

The parameters contained in the two-light-meson potentials such as mji, gab,R, Cab,R, 
^a^b',abi and by are determined by fitting experimental data. A particular choice of the 
potential, such as the number of R and contact interactions included, will be specihed later 
for each partial wave and for each of tttt and Kn interactions. 


B. Three-light-meson scattering model 


We now consider a case where three light-mesons are scattering. First, let us assume 
that the three mesons interact with each other only through the two-meson interactions 
discussed in the previous subsection. Because our two-meson interaction potential is given 
in a separable form, we can cast the Faddeev equation into a two-body like scattering 
equation (the so-called Alt-Grassberger-Sandhas (AGS) equation [22]) for a cTZ —)■ c'TV 
scattering. Here, TZ stands for either R or Tab, and Tab is a spurious “state” that is supposed 
to live within a contact interaction in a very short time, and decays (going to the left in 
equations) into the two light-mesons, ab. This degree of freedom is introduced merely for 
extending the AGS-type cR —)■ c'R' scattering equation used in Ref. [13] to include the 
contact two-meson interactions. Thus, the scattering equation for a partial-wave amplitude, 
{cn)i^P 'given as (see Fig. 2 for a diagrammatic representation) 


rjilJPT 


{p',p-,E) 


yC,JPT 


{p',P]E) 



yc",JPT 


ip', q; E) 


xGc"n’",c"n"iq, ’ 


(17) 


where JPT are the total angular momentum, parity, and the total isospin of the clZ system 
and they are conserved in the scattering. The cTZ state with the relative orbital angular 
momentum I is denoted by {cTZ)i] the allowed range for I is determined by JR and the 
spin-parity of TZ. The magnitude of the incoming (outgoing) relative momentum of the 
cTZ [c'TZ') state is denoted by p (p'). The driving term of the scattering is a partial-wave 
form of the so-called Z-diagram, E). The Z-diagram is a process in which 




FIG. 2. Diagrammatic representation of scattering equation FIG. 3. Z-diagram for cTZ —)■ 
of Eq. (17) for cR. dR! scattering. The gray blob represents dR! process via c-exchange. 
the dressed R, Green function in Eqs. (18)-(21). 



7^ —>■ c'c decay is followed by cc —?■ T^'-formation, as illustrated in Fig. 3. For a more 
explicit definition as well as the partial-wave expansion of the Z-diagram, we refer the 
readers to Appendix C of Ref. [13]; in particular, Eqs. (C10)-(C12) of Ref. [13] give an 
explicit expression for E) (note TZ=R and 7l'=R' in Ref. [13]) in which 

fdLn and dehned in Eq. (10) can be directly inserted. When 71 = r^'c {71' = the 
corresponding expression for the Z-diagram can be practically obtained by replacing 

in with w^,i dehned in Eq. (4). The Z-diagrams 

are known to have the moon-shape singularity [23] that prevents us from solving Eq. (17) 
with the standard subtraction method. Here we employ the spline method (see Ref. [23] for 
detailed explanations) to obtain numerical solutions from Eq. (17). 

In Eq. (17), we have also used the Green function, Gcn',cn{<ii E), for TZ and 7Z' which can 
be coupled through 7Z ah ^ 7Z'. It is given by 


[G-\q, =[E- E,{q) - E«(g)]<5«^« - (g, E - E,{q)) , 


for (7^,7^') = {R,R'), and 
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= -4/(g, E - E,{q)) for (7^, 7^') = (r„,, R') , 


(18) 

(19) 

( 20 ) 




rab,r^,y 


' — (T, 


LI 


,rab 


(g, E - Ec{q)) for (7^, 7^') = (r^b, r'^,y) , (21) 


where we have introduced the self-energies j^{q, E), and also quantities ■^^{q, E), 
which is either dimensionless or dimension of square-root of the energy, as dehned by 


'^^',r{Pi E) — 


ab 


mR/rriR 
Er'{p)Er{p) Jq 


°° 2 , 
q dq 


fR',ab{^)fab,R{^) 


^r''R{P^E) = J2 


ab 


rriR 

Er{p) 


q^dq 


^kirJP.E) = 


rriR! 

Er'{p) 


q^dq 


\/Mlb{q) +p‘^E- ^/M^^{q) + p‘^ + ie 
Mabjq) h^'i',abWab{(l)fL^,R{^) 

^Ml{q) +p^E- ^Ml{q) + p^ + ze' 
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-.( 22 ) 

(23) 

(24) 
(26) 
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FIG. 4. Vector iV) and pseudoscalar (P) 
mesons interaction potentials based on the 
hidden local symmetry model [21]. Vector 
mesons iVex) are exchanged. 


FIG. 5. Three-meson force not considered in 
this work, (a) Vector-meson exchange between 
a pseudoscalar-meson pair {ab) in s-wave and 
a pseudoscalar meson (c). (b) Rc interaction 
via meson-resonance (M*) excitation. 


with Mab{q) = Ea{q) + Eb{q), and J2ab two-meson states from R ^ ah decays. 

The kinematical factors in the expressions are from the Lorentz transformation to boost the 
7?.-at-rest frame to the cTZ center-of-mass frame. 

So far, we have considered the three-meson scattering due to multiple iterations of the 
two-meson interactions. Given the two-meson interactions, this is a necessary consequence 
of the three-body unitarity. In a three-meson system, however, there may be a room for 
a new mechanism that is absent in a two-meson system to play a role. We will refer to 
such mechanisms as a “three-meson force” hereafter. Diagrams shown in Fig. 4 can work 
as a three-meson force. These are interactions between a vector-meson and a pseudoscalar 
meson via a vector-meson exchange; for our particular application to —)■ , they 

are bare p-K, and bare K*-7r interactions. These mechanisms are based on the hidden local 
symmetry (HLS) model [21] in which vector and pseudoscalar mesons are implemented in a 
Lagrangian that has a symmetry under nonlinear chiral transformations. Expressions for the 
Lagrangian and the resulting interaction potentials of Fig. 4 are presented in Appendix A. 
These mechanisms in Fig. 4 along with the Z-diagram in Fig. 3 have been studied by Jansen 
et al. [24-26] to examine the vr-p correlation and its relevance to a soft ttNN form factor in 
a NN potential. There are also other possible mechanisms that can work as a three-meson 
force. We show some diagrams in Fig. 5 as examples. The diagram in Fig. 5(a) describes 
an interaction between a pseudoscalar-meson-pair (ah) in s-wave and another pseudoscalar 
meson (c) via a vector-meson exchange; this mechanism is also from the HLS Lagrangian. 
Meanwhile, in the diagram of Fig. 5(b), an R interacts with a pseudoscalar meson to form a 
resonance {M*), which is followed by a decay into an R' and a pseudoscalar meson. This is 
a familiar mechanism and often assumed in partial-wave analyses for meson spectroscopy. 

In this work, we consider the vector-pseudoscalar interactions shown in Fig. 4 in our 
analysis of the —)■ K~7r'^7r'^ decay to study their relevance. Thus, the scattering equation 
in Eq. (17) is modified by adding the new mechanisms to the Z-diagrams: 


P'^ ^ ip'^ p; p) in Eq. (17) , (26) 

where the added term, (p', p), is in the partial-wave form for which we give explicit 


10 



















FIG. 6. Diagrammatic representation of D'*'-decay into three pseudoscalar mesons a, b, c in 
our coupled-channel model, (a) The isobar-type diagram, (b) The rescattering diagram. The 
amplitude T' is from the scattering equation represented by Fig. 2. The gray blob represents the 
dressed IZ Green function. 


expressions in Eqs.(A8), (A12), and (A17). In Eq. (26), R and R' are the lightest spin-1 bare 
states of either (/, S'[strangeness]) = (1,0) or (1,5') = (1/2,—1) which we denote “p” and 
“i?*”, respectively. For our particular application to the iT“7r+7r+ decay, we include 

{VP,V'P',Ve^)= (V’iT,‘y’J?,p), (‘y/ir,“iT*”7r,A:*), (“ir*”7r,“iT*”7r,p) for the_diagram 
Fig. 4(a), and {VP,V'P',Ve^)= K,K*), (‘y’A',“iT*”7r,a;), (“A:*”7r,“A:*”7r,A:*) 

for the diagram Fig. 4(b). On the other hand, we leave examination of mechanisms such 
as those shown in Fig. 5 to future work for the following reasons. As we emphasized in the 
introduction, even effects of multiple scattering due to the two-meson force on the H-decay 
has still not been studied in a realistic setting. In this situation, studying the relevance of 
the three-meson force to the H-decay is indeed in an exploratory level, and thus a reasonable 
starting point would be to include it in the most important channel. As we will see, the 
vector-pseudoscalar {p-K) channel plays a very important role in the rescattering process, 
and therefore it would be good to study the new mechanisms of Fig. 4 in this channel at 
hrst. Also, regarding the mechanism in Fig. 5(b), no relevant meson-resonance of spin-0, 
1=3/2, S=—l is known in the D-meson mass region, and thus we do not need to include it 
for the moment. 


C. D+ — K 7r+7r+ decay amplitude 

In our coupled-channel model, the decay amplitude for —)■ is given by 


cyclic 

^iC“7r+7r+,D+5 PttJ*" 5 PttJ 5 ^ ^ ^{afo)c,D+(pa? Pc;-S') 5 (2^) 

(abc) 


where we have introduced the cyclic summation that takes the sum over abc = K tTiTT^, 
111 712 K~, 

T[ah)c,D+{.Pa,Pb,Pc]E) = ^ Pb)Gc7e,c7e'(Pc ^)f (Pc, ^), (28) 

71 , 71 ', 

where is the 2 ;-component of the spin of TZ, and Gc 7 i,c 7 Z' is the Green function that 
has been dehned in Eqs. (18)-(21). This decay amplitude in Eq. (28) is diagrammatically 
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s ^ 

represented in Fig. 6. The symbol denotes a 7^ —> a6 decay vertex function which is 
explicitly given as 


fL^niPa, Pb) 

fab,APa,Pb) 


mREa{q)Eb{q) 

Eji{pji)EaiPa)Eb{pb) 


{tafjbttll, tl + tl)YLp4q)fabM for 7^ = R{29) 


Ea{q)Eb{q) 

Ea ip a) Eb {pb j 


{tatltbtl\I, tl + tl)YLp4q)wj^b (q) for n = rab , (30) 


where f^bR ^rid have been dehned in Eqs. (9) and (4), respectively; (jimij 2 "i 2 |TM) is 
the Clebsch-Gordan coefficient, and ta is the isospin of meson a and is its .2 component. 
The kinematical factors in the equations are from the Lorentz transformation to boost the 

a6-pair center-of-mass (CM) frame to the total CM frame. The momentum q is the relative 

_ 

momentum of the a6-pair in their CM frame. The dressed TZ'c vertex, T^^', , has 

also been introduced in Eq. (28), and it is explicitly given by 


r" 


J^,^rAPc,E) = t^+tl tR:\T tl+tl+t:) 

1,1^ 


PT 

xY, 


LP 




{Pc,E), 


(31) 


where Sd{= 0) is the H-meson spin, and = 3/2 for —)■ K~7r~^7r~^. We sum over 

the parity (P) and total isospin (T) of the hnal hadronic states because the weak P-decay 
does not conserve them. The last factor in the above equation is given by 

/ CX) 

dpc pI (,,7^,,)^, (p„ Per, E) 

xGc'n'\c''R."'iPc', E)E^f^,„^^^ R+{pc') , (32) 

where is the partial-wave amplitude for c'TZ” —?■ cTZ' scattering obtained by 

solving the coupled-channel scattering equation, Eq. (17). The hrst term on the rhs corre¬ 
sponds to the isobar-type contribution [Fig. 6(a)] while the second term is the contribution 
from the rescattering [Fig. 6(b)]. The quantity E^J^^ ^+{p) is the bare (Ec)i vertex 

function for which we choose a parametrization. 




An C(cnh 






2 \ 2+{1/2) ^ ^ 

P 


(2Kfr V?' + (HYP 


— , (33) 


with 


— 


rriR 


‘2Er{p) 


for 77. = P , An = 1 for 77. = Tab ■ 


(34) 


I„ Eq. (33), C"l, pp , and are the coupling, phase, and cutoff, respectively, and 

they will be determined by htting Dalitz plot distribution data. The couplings are 

nonzero only when IFo —st^I < I < Sn + sn- The parametrization used in this work [Eq. (33)] 
is different from the one used in Ref. [13] in choosing the kinematical factor. 
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The isobar-model-type amplitude [Fig. 6(a)] for the —)■ K 7r+7r+ decay, , is 

obtained from the above equations (27)-(34) by just dropping the second term of the rhs of 
Eq. (32). 

The procedure to calculate the Dalitz plot distribution from the decay amplitude of 
Eq. (27) is explained in detail in Appendix B of Ref. [13], and we do not repeat it here. 


III. ANALYSIS AND RESULTS 

Now we apply the coupled-channel formalism discussed in the previous section to analyses 
of data. First we determine parameters in the two-pseudoscalar-meson scattering model by 
analyzing experimental data for tttt and nK scatterings. Then we extract resonance parame¬ 
ters from amplitudes of the two-meson interaction model. This two-meson interaction model 
is a basic ingredient for the three-meson scattering model. In the subsequent subsection, we 
analyze the —)■ decay in a realistic setting. 


A. Two-pseudoscalar-meson scattering 

For studying the —)■ decay in our coupled-channel framework, the tttt and 

ttK scatterings of E < 2 GeV are relevant. We will determine the model parameters of our 
TTTT and ttA scattering models, i.e., rriR, gab,R, and Cab,R in Eqs. (l)-(3), (8), and 

(10), by fitting empirical scattering amplitudes for E < 2 GeV. 


1. ttK scattering 

We analyze the tiK scattering amplitudes from the LASS experiment [6, 9]. For our 
application to the —)■ decay in the next section, we determine the model pa¬ 

rameters for {T,/} = {0,1/2}, {0,3/2}, {1,1/2}, {2,1/2} partial waves. We explain details 
of our t[K scattering model for each partial wave. For the {L, /}={0,l/2} wave, we consider 
TrK-rj'K coupled channels because the r]'K channel is known to play a significant role while 
rjK does not in this partial wave. We include two bare R states supplemented by a contact 
ttK^tiK interaction. For the {L,/}={0,3/2} wave, we consider a contact inter¬ 

action only. For the {L,/}={!,1/2} and {2,1/2} waves, we consider coupling of ttK and 
effective inelastic channels; masses of the two “particles” in the inelastic channel, denoted by 
and , are also htted to the data. We include three bare R states for {L, /}={!,1/2} 
while a single bare R state for {2,1/2}. We present the ttK model parameters determined 
by the hts in Table VII of Appendix B. 

We present the quality of the hts to the empirical partial-wave amplitude [6] of vr+JF” 
L=0 partial wave and {L,/}={!,1/2},{2,1/2} partial waves in Fig. 7 where phases (upper 
panels) and modulus (lower panels) of the amplitudes are shown. Also the elastic scattering 
phase shifts for the {L,/}={0,3/2} partial wave calculated with our model are compared 
with the data [9] in Fig. 8(left). The values obtained in the hts are tabulated in Table I. 
The s-wave vr+iF” amplitude is calculated by linearly combining the {L, /}={0,l/2}, {0,3/2} 
partial-wave amplitudes. Overall, as seen in the hgures, we obtain a reasonable description 
of the data included in the ht. However, one notices that the model has a sudden change and 
deviation from the data in the phase for the {L,/}={!,1/2} partial wave at E ~ 1.3 GeV. 
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600 1000 1400 1800 

E (MeV) 



FIG. 7. (Color online) Phase (upper) and modulus (lower) of the amplitudes for the ttK scat¬ 
terings: (Left) L=Qioi'K~^K~ ; (Center) {L,/}={!,1/2}; (Right) {L,/}={2,l/2}. Data are taken 
from Ref. [6]. 




E (MeV) E (MeV) 

FIG. 8. (Color online) Phase shifts {6li) of the ttK and tt-tt scatterings: (Left) ttK phase shifts 
for {L,/}={0,3/2} partial wave; data are from Ref. [9]. (Right) Trvr phase shifts for {L,/}={0,2} 
partial wave; data are from Refs. [27, 28]. 


This is perhaps an artifact of our model that has the threshold where the effective inelastic 
channel opens. Fortunately, the magnitude of the amplitude is rather small around this 
energy so that the deviation in the phase will not give a significant impact on observables 
calculated with this model. 

From the ttK amplitudes of the model determined above, we extract resonance pole 
positions by the analytic continuation [31, 32] as shown in Table II. We present poles below 
Re[i?] < 2 GeV and ]Im[Fl]] < 0.25 GeV. We can consistently identify most of the extracted 
poles with the corresponding particles listed by the Particle Data Group (PDG) [30] as shown 
in the table. For the {L, J}={0,l/2} partial wave, our model has a pole at 0.71 — 0.23f GeV 
that corresponds to the so-called k whose mass is 682 ± 29 MeV, and width 547 ± 24 MeV 
in the PDG listing. Also we find two poles at Re[i?] ~1.4 GeV on different Riemann sheets. 
These two poles are associated with a single resonance [iFo(1430)] that is split by coupling 
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TABLE I. values for the fits to empirical -kK and tttt partial-wave amplitudes. Each entry in 
the row indicated by of parameters” is the number of model parameters adjusted to fit the 
data. 


{L,I} 

{0,1/2} 

{0,3/2} 

ttK 

{1,1/2} 

{2,1/2} 

{0,0} 

{0,2} 

TTTT 

{1,1} 

{2,0} 


304 

158 

344 

183 

274 

114 

291 

276 

# of data points 

84 

19 

84 

64 

148 

42 

130 

130 

^ of parameters 

12 

2 

17 

7 

12 

3 

10 

5 

xVd-o.f 

4.2 

9.3 

5.1 

3.2 

2.0 

2.9 

2.4 

2.2 


TABLE 11. Pole positions of the ttK partial-wave amplitudes in the complex-energy plane. A 
partial wave is specified by the orbital angular momentum (L) and the isospin (I). Poles below 
Re[E] < 2 GeV and |Im[E]| < 0.25 GeV are presented. Roman numerals in the square brackets 
specify the Riemann sheet on which each of the poles exist. We use the convention defined in, 
e.g., Ref. [29], to specify each of the Riemann sheets, I-IV. Each of the states is identified with the 
corresponding particle name used in the PDG listings [30]. 


L I 

0 1/2 0.71 

1 1/2 0.90 

2 1/2 1.42 


Pole positions (GeV), [Riemann sheet]. Name 


0.23i [II] K 

0.025i [II] A* (892) 
0.055i [III] A|(1430) 


L44-0.14i [II] A*(1430) 
1.46 - 0.25i [III] Ad (1430) 
1.28 -0.058i [III] A* (1410) 


1.88-0.13Z [III] Ad (1950) 
1.66 -0.088i [III] A* (1680) 


to rj'K channel. For the {L,/}={!,1/2} partial wave, our model has the well-established 
A*(892). Also in the same partial wave, there is a pole at 1.28 — 0.058z GeV that is a bit 
off the A* (1410) resonance parameters from the PDG. 


2. TTTr scattering 

We perform an analysis of tttt scattering data with our coupled-channel model in a way 
similar to the analysis of ttK data in the previous section. Although we only need a tttt model 
for the {L,/}={!,1},{0,2} partial waves for our coupled-channel analysis of the —)■ 
A“7r+7r+ decay, we present here our tttt model for all major partial waves for a future 
reference. We consider tttt-AA coupled channels for all partial waves except for {L, /}={0,2} 
where only the elastic tttt channel is taken into account. Regarding details of our model 
for each partial wave, we include two bare R states supplemented by a contact tttt—j-tttt 
interaction for the {L,/}={0,0} wave. For the {L,/}={!,1} and {2,0} waves, we include 
two bare R states and a single bare R state, respectively. Finally for the {L,/}={0,2} 
wave, we consider a contact tttt— )-7r7r interaction only. We present the tttt interaction model 
parameters determined by the fits in Table VIII of Appendix B. 

We present the quality of the hts in Figs. 9 and 8(right) where phase shifts and inelastic¬ 
ities are shown. The values obtained in the fits are tabulated in Table I. As seen in the 
figures, we obtain reasonable fits to the data from Refs. [18-20]. Although the quality of 
the fits is not much improved from those we obtained in Ref. [13], the purpose of updating 
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FIG. 9. (Color online) Phase shifts (upper) and inelasticities (lower) for the tt-tt scattering: (Left) 
{L,/}={0,0}. (Center) {L,/}={!,!}• (Right) {L,/}={2,0}. Data are taken from Ref. [18-20]. 


TABLE III. Pole positions of the vTTr partial-wave amplitndes in the complex-energy plane. The 
other features are the same as those in Table 11. 


L 

I 



Pole positions (GeV), 

[Riemann sheet 

, Name 

0 

0 

0.41 - 0.24i 

[11] 

a 

1.02 - O.OOH 

[11] 

/o(980) 

1.38 - 0.15i 

[Ill] 

/o(1370) 

1 

1 

0.77 - 0.079i 

[11] 

p(770) 

1.01 - 0.083i 

[III] 

p(1450) 

1.58 - 0.16i 

[III] 

p(1700) 

2 

0 

1.25 - 0.099i 

[III] 

/2(1270) 

— 



— 




the TTTT scattering model is to take account of the proper kinematical factor, as mentioned 
below Eq. (10). The nonzero inelasticities from our model are due to the coupling to the 
KK channel of our coupled-channel model. We note that the KK channel in our model 
effectively simulates all inelastic channels in which the true KK channel is a major compo¬ 
nent, because we did not include tttt KK and KK ^ KK data in our analysis, and also 
because we did not include the dvr channel in the model. From the determined tttt partial- 
wave amplitudes of our model, we extract resonance poles as presented in Table III. Most 
of the extracted poles are consistently identihed with the counterparts in the PDG listings, 
as shown in Table III. A difference from the PDG value is found for the width of /o(980); 
our model has a rather small width (~2 MeV) in comparison with the PDG average (40- 
100 MeV). This difference was also found in the model used in Ref. [13], and possible sources 
of the difference were discussed there. Also the second resonance in the {L, /}={!,1} partial 
wave perhaps does not correspond to p(1450). However, an effect of this resonance pole on 
the amplitude seems to be very small, and our model reproduces the empirical amplitude 
very well. These differences from the PDG listings in the pole positions could be due to the 
fact that we employ the simple parametrization for the two-pseudoscalar-meson interactions 
in order to make our three-meson scattering model relatively easy to handle. The behavior 
of the amplitudes on the unphysical (complex energy) region may be different from those 
of more sophisticated two-meson interaction models cited in the PDG. Thus, we would not 
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mean to claim that the pole positions presented on Tables II and III have a comparable 
reliability to those from the more sophisticated analyses. The tables are just for showing 
the properties of the amplitudes used in our analysis of the —)■ decay, and thus 

we do not quote errors for the pole positions. Yet, the tttt and •nK amplitudes our model 
generates are reasonable on the real physical energy axis, and should be good enough for 
our purpose, that is, a coupled-channel analysis of the —)■ ii'“7r+7r+ decay with realistic 

TTTT and ttK amplitudes. 


B. Analysis of —)■ K Dalitz plot 

Now we will perform a partial-wave analysis of pseudodata for Dalitz 

plot distribution using our coupled-channel model. In what follows, we explain setups of our 
models used in the analysis. Then we discuss how we prepare pseudodata and our analysis 
procedure, which is followed by numerical results. 


1. Model setup 


In our coupled-channel framework, D’''-meson decays into TZc channels, followed by mul¬ 
tiple scatterings due to the hadronic dynamics, leading to the hnal state. This 

process is expressed by Eqs. (28), (31), and (32); with the symmetrization, the decay am¬ 
plitude is given by Eq. (27). We consider the following 11 TZc coupled channels in our full 
calculation: 


{TZc} = {Rfn, i?2^7r, TT, Rfw, Rl‘^K, RfK, 


TT 

’ ' TTTT 


/?} , (35) 


L 2/ 

where R^ ’ stands for the i-th bare R state with the spin L and the isospin /; when I is an 
integer (half-integer), it is understood that this R state has the strangeness S'=0 {S=—l) in 
this paper. Thus, R^^, Rj^, Rf^, are seeds of Kq, K*, K 2 , p resonances, respectively. In 
our model, these resonances are included as poles in the unitary scattering amplitudes, 
is a “state” associated with a contact interaction in a partial wave of L and /, as has been 
introduced in Sec. IIB. Most of the partial waves associated with these channels have been 
considered in the previous Dalitz plot analyses of the D+ —)■ A“7r+7r+ decay [2, 3, 5, 8,10, 11]. 
However, the {L, J}={0,2} partial wave associated with the channel was considered 

only in the CLEO analysis [11]. Also, the {L,/}={0,3/2} partial wave associated with the 
r^^TT channel was explicitly considered only in the FOCUS analysis [8], but other MIPWA 
can implicitly take account of this partial wave. The {L,/}={!,1} partial wave associated 
with the Rj^K channel has not been considered in the previous analyses. This channel can 
contribute to —)■ K~7r~^7r~^ only through the coupled-channel dynamics, and therefore 

it does not show up in isobar-type models. A possible important role of this channel was 
hinted by the Brazilian group [15, 17], as stated in the introduction. We note that, unlike 
most of the previous isobar-type analyses, we do not include a flat interfering background 
amplitude. With the coupled channels considered in this work [Eq. (35)], the hnal hadronic 
system has the total spin J=0, parity P=-|-l, total isospin T=3/2, and 1=sti. We ht 
Dalitz plot pseudodata for by adjusting parameters associated with the 

TZc vertex function in Eq. (33). Among the parameters, and we 

£x GeV for all D~^ — {cTZ)i vertices, and set Then, we ht 


17 


the pseudodata by adjusting the other and under a requirement that all , 

which belong to the same partial wave characterized by {L, /}, share the same phase 
With this requirement, when the hadronic rescattering is turned off, the Watson theorem 
is satished, up to a slight violation due to the fact that our model takes account of the 
center-of-mass motion of the scattering two mesons. More specifically, the pc-dependence of 
the phase of Genin' iu Eq. (28) leads to the slight violation. 

The hadronic rescattering processes are described by the scattering amplitude, (c 7 i)i ’ 

dehned in Eq. (17), and also by the TZe Green function dehned in Eqs. (18)-(21). The main 
driving force of the rescattering processes is the two-pseudoscalar-meson interactions that 
have been hxed in the previous section by fitting the empirical ttK and tttt scattering am¬ 
plitudes. The two-pseudoscalar-meson interactions enter the scattering equation [Eq. (17)] 
as the Z-diagrams and also through the R self-energies and in Eqs. (22)-(25). Given 
the coupled channels specihed above and the two-meson interactions, we can find all the 
contributing Z-diagrams that induce channel couplings. In our analysis, we consider all of 
the contributing Z-diagrams that contain tittK and KKK intermediate states. In addition, 
we include the three-meson force based on the HLS model and, as mentioned below Eq. (26), 
we have totally six diagrams of this type. As described in Appendix A, once two coupling 
constants are fixed by the p — )■ tttt and cu —)■ vrp decay widths, all the other couplings are 
fixed by SU(3) and the HLS. We use the same form factor [Eq. (A15)] for all the different 
diagrams of the three-meson force. The cutoff, Asmf, in the form factor is determined by 
fitting the pseudodata. 

In our analysis of —)■ Dalitz plot pseudodata, we will basically use three 

models. The hrst one is the “Full model” that contains all the dynamical contents described 
above. The second one is the “Z model” for which the rescattering mechanism is solely due to 
multiple iteration of the two-pseudoscalar-meson interactions in the form of the Z diagrams 
and TZ propagators. Thus the Full model is obtained by adding the three-meson force to 
the Z model. The third model is the “Isobar model” that does not explicitly contain the 
rescattering process. The decay amplitude for the Isobar model has been described at the 
end of Sec. IIG. This Isobar model is still different from most of the isobar models used in the 
previous Dalitz plot analyses of —)■ in that all two-pseudoscalar partial-wave 

amplitudes are unitary, and ht well the empirical amplitudes in the relevant energy region; 
the Watson theorem is also maintained in the sense discussed in the previous paragraph. 
Finally we remark that the two-pseudoscalar-meson interactions, that have been determined 
in the previous sections, will not be adjusted to ht the D+-decay pseudodata. This is in 
contrast with most of the previous analyses where some of Breit-Wigner parameters were 
also adjusted along with D+ —)■ Rc vertices. 


2. Data and analysis method 

We create reasonably realistic pseudodata of the —)■ Dalitz plot from the 

isobar model of the E791 Gollaboration [5]. The E791 Gollaboration obtained the isobar 
model through their partial-wave analysis of the —)■ K~7r~^7r~^ Dalitz plot of 15,079 

events, among which 94.4% were determined to be signals. In generating pseudodata, we 
take a procedure similar to Sec. V of Ref. [33] where Dedonder et al. created pseudodata 
for —)■ Kg7r^7r~ using the isobar model of the BABAR Gollaboration. We start with 

a grid 400x400 squared cells covering all kinematical region of the Dalitz plot distribution 
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FIG. 10. (Color online) The Dalitz plot distribution for D~^ —)• K~'k^'k^ in 40x40 cells. The 
left panel shows our pseudodata generated with the isobar model of the E791 Collaboration [5]. 
The right panel is the counterpart from the Full model that has been fitted to the pseudodata. An 
explanation for how the pseudodata are generated is found in Sec. IIIB2. 


with as x-axis and another as |/-axis; denotes the squared invariant 

mass of the K~ti^ pair. Each cell is given by the E791 isobar model a value of the Dalitz 
plot distribution at the center of the cell. Then, the values of the Dalitz plot distribution in 
10x10 adjacent cells are summed to obtain 40x40 cells, each of which has the value of the 
partially integrated Dalitz plot distribution. The width of each cell is 0.0649 GeV^. In the 
E791 analysis [5], 40x40 cells were used to perform their MIPWA and thus are a reasonable 
size also in our analysis. The Dalitz plot distribution value in each cell is multiplied by 
a common normalization constant and then is round off to be an integer; the common 
normalization constant is chosen so that the sum of the round-off values of all the cells 
coincides with 15,079 x 94.4% ~14,234 , the number of signals of the E791 experiment. In 
this way, we have generated pseudodata for the D+ —)■ iT“7r+7r+ Dalitz plot, as presented 
in the left panel of Fig. 10. 

Next task is to analyze the above pseudodata with our model. Again, Ref. [33] serves as 
a useful reference to hnd an analysis procedure. We calculate the Dalitz plot distribution 
using the decay amplitude of Eqs. (27) and the formulas given in Appendix B of Ref. [13]. 
In each cell of the 40x40 grid, we integrate the Dalitz plot distribution from our model; the 
overall normalization is chosen so that the integral over all the kinematical region of the 
Dalitz plot distribution becomes equal to 14,234, the number of signals for the E791 data. 
In this way, we have the number of events (a real number) in each of the cells, and can 
compare it with the counterpart in the pseudodata. If a given cell of the pseudodata has 
the number of events less than 5, then the cell is merged with the adjacent cell in the same 
x-axis to be a larger cell. This grouping is repeated until the cell contains more than or 
equal to 5 events. The same grouping is also applied to the event samples from our model. 
Thus, the total number of effective cells, each of which contains more than or equal to 5 
events of the pseudodata, is 725. 


19 








TABLE IV. The total x^-values and x^/d.o.f. from the Full, Z and Isobar models obtained by fitting 
the —)• K~T:~^Tr~^ Dalitz plot pseudodata. The Z model without couplings to is 

labeled by “Z (without /?)”. 



Full 

Z 

Isobar 

Z (without p) 


157. 

119. 

303. 

216. 

xVd-o.f. 

0.22 

0.17 

0.42 

0.30 


The quality of the £t can be quantihed by calculating yY defined by 


A 


2 




J^th _ 

J J 


2 




(36) 


where j labels each cell, and is the number of events in the j-th cell from our 

model (pseudodata). The error of the pseudodata in each cell is assigned by = 



We will perform the least x^”fit fo fhe pseudodata with the CERN program library 
MINUIT. Errors for the htting parameters are estimated in a standard manner as follows. 
First we calculate Hessian matrix, Hij, dehned by 




2 dOidOj 


{e}={e} ’ 


(37) 


where x^ is given in Eq. (36); 9i is one of the htting parameters and {6} is a set of the 
htting parameters at the minimum of x^- Then the error matrix is given by the inverse of 
the Hessian, Eij = the error for the parameter 6^ is assigned by y/Eu- An error 

6X for a quantity X such as a ht fraction is estimated by the error propagation formula: 


[M]"=i; 

id 






{e}={e} 




{ 0 }={ 0 } 


(38) 


3. Numerical results and discussions 

We performed the least x^"fif to fho pseudodata following the procedure explained in the 
previous subsection. We used the three models; the Full, Z, and Isobar models. All the 
parameters and their statistical errors determined by the hts are tabulated for each of the 
models in Table IX. The Dalitz plot distribution from the Full model is shown in the right 
panel of Fig. 10. Comparing with the left panel of the pseudodata, a diherence is hardly 
discernible to the eye. The situation is the same for the Z and Isobar models. The quality 
of the ht can be quantihed by calculating the x^-values of Eq. (36) that are presented in the 
second row of Table IV. The x^-value of the Z model is the smallest, and the Full model 
comes in second. These models are signihcantly better than the Isobar model in the ht 
quality. In the third row of the table, we also show x^/d-o.f. to assess if the better x^ is 
simply due to more degrees of freedom in the hts. The number of the htting parameters for 
the Full, Z, and Isobar models are 16, 15, and 12, respectively, as can be found in Table IX. 
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FIG. 11. (Color online) The Xj distributions. The left, middle, right panels are from the Full, Z, 
and Isobar models, respectively. The bins with Xj ^ ^ s-d given the same color (red); the same 
applies to the bins with xj < 0.5 (dark blue). 




FIG. 12. (Color online) The K 7r~^ (7r'*'7r'*“) squared invariant mass spectrum of the Full model 
compared with the pseudodata from the E791 isobar model in the left (right) panel. 


The number of bins at which x^ is calculated is 725. Thus we obtain x^/d.o.f. as shown in 
Table IV. Therefore, the ranking of the £t quality is still in the same order as far as x^ is 
concerned. 

In order to see the quality of the hts more clearly, we show the Xj distribution at all 
of the bins in Fig. 11. These hgures show that all of our models ht the pseudodata rather 
precisely. Most of the bins are htted with Xj < 0.5, and Xj exceeds 1 at only a small number 
of the bins. Yet again, the Full and Z models clearly show a better performance in the £t 
than the Isobar model does. The quality of the £t can also be shown in the projection of the 
Dalitz plot distribution onto the or distributions, as presented in Figs. 12 for 

the Full model as a representative. 

What is the main reason that the Full and Z models £t the pseudodata better than the 
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Isobar model ? Because only the former models have the Rj^K (i.e., p'^K^) channels and 
the hadronic FSI, one may guess either or both of these dynamical contents are responsible. 
To address this point, we introduce the “Z (without p)” model, that is, the Z-model with 
couplings to the R]^K channel turned off. By htting the pseudodata with the Z (without p) 
model, we examine if the rescattering effects can improve the The value is shown in 
Table IV, and is signihcantly better than that of the Isobar model. We note that the Isobar 
model and the Z(without p) model have the same number of adjustable parameters in the 
hts. Thus, the quality of the £t is improved just by including the rescattering due to the 
Z-diagrams, and the inclusion of the Rj^K channels further improves the fit substantially. 

We remark that we obtained the reasonable hts to the pseudodata without adjusting the 
parameters associated with the two-pseudoscalar-meson interactions. On the other hand, as 
mentioned already, most of the previous analyses varied some Breit-Wigner parameters in 
htting their D~^ —)■ K~7r~^7r~^ Dalitz plot data. A common claim in those analyses [2, 5, 8 , 11] 
was that the width of A’q(1430) obtained in their hts, ~170 MeV, was signihcantly smaller 
than those from the PDG and the LASS analysis, ~270 MeV. In our model, there are two 
poles associated with A^^(1430) as shown in Table II, and they are on diherent Riemann 
sheets, the branch point of which is the p'K threshold. This two-pole structure is a coupled- 
channel ehect. One of the poles is close to the PDG value, while the other one has a rather 
broad width. Thus, in our analysis, we did not need JPq(1430) with a narrower width to 
obtain the reasonable hts to the —)■ Dalitz plot pseudodata. 

Next we present partial-wave decay amplitudes dehned by Eq. (28) at g • Pc = 1 {.Q'- 
relative momentum of the a 6 -pair in their GM frame) and the summation 'Ylnw being 
replaced by Yhn' where {T,/} specify the partial wave. These ampli¬ 

tudes correspond to Eq. (3) of Ref. [5] at cos 6 ' = 1, and in the same reference, the E791 
Gollaboration presented their amplitudes. Thus we can compare our amplitudes with those 
from the E791 MIPWA. We denote a partial wave by “( 06 )^ 0 ” in which a two-pseudoscalar- 
meson pair ab has the total angular momentum L and the total isospin I. When ab is written 
with a charge state, (afe)£ is understood to be the total isospin I state projected onto the 
particular charge state, e.g., (iP“ 7 r+)g“^'^^ 7 r+. We note that a partial-wave amplitude of 
a coupled-channel model generally contains all the partial waves considered in the model 
as intermediate states. Therefore, we refer to an amplitude with (a&)^c in the hnal state 
as the {aby^c amplitude. The partial-wave amplitudes for (iP“ 7 r''')^~^'^^ 7 r+, (iP“ 7 r’'')p“^'^^ 7 r’'', 
(A“ 7 r’'')jj“^'^^ 7 r+, (iP“ 7 r+)g“^'^^ 7 r+, and K~ are presented in Figs. 13(a)-(e), respec¬ 

tively. The partial-wave amplitudes from the Full, Z, and Isobar models are shown with their 
error bands as a function of Mab = {Ea -l- EbY — {pa + PbY- As a whole, the Full and Z 
models are similar in the amplitudes while the Isobar model is rather diherent, particularly 
in (A'“ 7 r''')^“^^^ 7 r''', , and {7r~^Tr~^)g"'^K~. Because the models maintain the 

Watson theorem when the hadronic rescattering is absent (see the note in Sec. IIIBl), the 
phases (the upper panels of Fig. 13) from the Isobar model in the elastic region are essen¬ 
tially the same as those in Figs. 7 and 8 up to overall constant shifts. The diherence in the 
Mx- 7 r+-dependence of the phases between the Isobar model and the Full (or Z) model is 
purely the ehect of the hadronic rescattering. To put it the other way around, the Isobar 
model does not have a degree of freedom to change the phases using the rescattering in order 
to obtain an optimal ht, which might have led to the rather diherent solution. 

In Fig. 13(a), we also show the {K~7r~^)s7r~^ amplitude from the E791 MIPWA [5], denoted 
by [(A'“ 7 r’'') 57 r''']®’jp\^^ hereafter. The modulus in the hgure is, in their notation, c x F^{y/s), 
and numerical values for c and are tabulated in Table III of Ref. [5]. Interestingly, 
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FIG. 13. (Color online) Phase (npper) and modulus (lower) of the partial-wave ampli¬ 
tudes for (a) , (b) (K“7r+)p'~^'^^7r+, (c) {K~, (d) {K~, (e) 

{'k'^tt'^Y^'^■ The red solid, the blue cross-hatched, and the green bordered bands are from the 
Full, Z, and Isobar models, respectively; the band widths represents the errors. For the dehnition 
of the amplitudes, see the text. In the panel (a), we also show the (iir“7r+)5'7r'’' amplitude from 
the MIPWA of the E791 Collaboration [5] as the black squares with error bars. 


the dependence of the phases from the Full and Z models are in a very good agree¬ 

ment with those of the MIPWA for Mk--k+ ^1-5 GeV. The Full model also agrees with the 
modulus of [(A'“7r+)57r+]®'(p\,y^ in the elastic region. Because [(A'“7r+)57r+]®'(p\,y^ implicitly 
contains all the partial-wave amplitudes other than (iF“7r+)p~^^^7r’'' and (iF“7r’'')^^'^^7r’'', 
for a more meaningful comparison, we might need to compare it with the coherent sum of 
the j (7r+7r+);5^^iF“ amplitudes of our models. How- 
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ever, because the (vr’''7r+);5 amplitude depends on M^+^+ while the others on Mk--k+i 
they cannot be simply summed to obtain the counterpart of Still, we 

conhrmed that the Mi^-^+-dependence of the phase of the amplitude of the 

Full and Z models does not significantly change even after the {K~7r~^)g~ ' tt’*' amplitude 
is coherently added. Thus, our coupled-channel models explain the gap between the phase 
of and that of the LASS {K~TT^)g~^^^TT^ amplitude in a way qualita¬ 

tively different from the previous explanation. Edera et ah [7] and the FOCUS iC-matrix 
model analysis filled the gap with a rather large (more than 100%) destructive interference 
between the and amplitudes, without an explicit consider¬ 

ation of the hadronic rescattering. Our coupled-channel models, on the other hand, hll the 
gap with the hadronic rescattering, and have a moderately destructive interference between 
the and amplitudes (see Table V). 

We notice in Fig. 13 that the and (7r’''7r’'');^^^iF“ partial-wave ampli¬ 

tudes have relatively large errors. Because we analyzed the data of the single charge 
state, it may be difficult to separate the different isospin states with a good precision. 
This situation would be improved by analyzing data of different charge states, i.e., D~^ 

and —)■ in a combined manner. We also note that even though the 

phases of the {K~Tr~^)jr ' 7r+ amplitude from the Isobar model have very large errors for 
Mk-tt+ ^1-5 GeV, this is simply because the absolute value of the amplitude is very small. 

Now we look into the fraction of each channel’s contribution (fit fraction). In an isobar 
model that describes a Zl''’-meson decay as D~^ —)■ —?■ abc where R is expressed by a 

Breit-Wigner function, the ht fractions of different Rc contributions can be straightforwardly 
dehned, and are often presented in the previous analyses. However, in a model like ours where 
the resonances are described as poles in unitary scattering amplitudes, the ht fraction of 
a certain resonance contribution is not so straightforwardly dehned, because the scattering 
amplitude generally contains more than a single resonance, as we have seen in Tables II 
and III. Furthermore, the amplitude also contains nonresonant contributions. There is no 
unambiguous way to single out a certain resonance contribution. Therefore, we present the 
ht fractions of contributions from diherent partial-wave amplitudes. We calculate the (a&)£c 
partial width using the partial-wave amplitudes presented in the above paragraph, with the 
angular [q ■ pc) dependence being restored and the Bose symmetrization being taken into 
account as in Eq. (27). The ht fraction is then naturally dehned as the {ab)\^c partial width 
divided by the total width. 

The ht fractions dehned in the above paragraph are presented in Table V for the Full, Z, 
and Isobar models. The incoherent sum of the ht fractions in the last row is not necessarily 
100 % because of the interferences between diherent contributions in the diherent rows of the 
table. We can see that the Full and Z models agree fairly well within the errors, although the 
errors are rather large for the and {7r~^7r~^)g"‘^K~ ht fractions. On the other 

hand, the Isobar model has quite diherent and (7r'''7r’'')^^^iF“ ht fractions 

that are rather large, leading to the large incoherent sum (367%). This indicates a very 
destructive interference between the and partial waves. These 

results are consistent with what we can expect from the partial-wave amplitudes shown 
in Fig. 13. We also show the coherent sum of the (iF“7r+)g“^'^^7r+, (JF“7r+)g“^^^7r+, and 
ht fractions, labeled “S-waves”. Interestingly, all three models have consistent 
“S-waves” ht fractions within the drastically reduced errors. Therefore, the data used in our 
analysis can constrain the “S-waves” ht fraction rather well, while they cannot well constrain 
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TABLE V. Fit fractions (%) from each of the models fitted to the —)• Ar“7r'*'7r'*‘ pseudodata. See 
the text for the definition of the fit fraction. The “S-waves” fit fraction is the coherent sum of the 

and fit fractions. The numbers in the parentheses 

are not included in the “Sum” in the last row. The fit fractions from the E791 Isobar model [5], 
the FOCUS A-matrix model analysis [8], and the CLEO QMIPWA [11] are also presented. The 
fit fractions in () are obtained by the incoherent sum of different resonance contributions in the 
same partial wave, and the errors are added in quadrature. 



Full 

Z 

Isobar 

E791 [5] 
isobar 

FOCUS [8] 

A-matrix 

CLEO [11] 
QMIPWA 

1 

+ 

lo 

+ 

95.9 ±7.3 

78.7 ±4.3 

68.9 ±4.4 

33.8 ± (10.8) 

207.25 ±25.45 

— 

1 

+ 

+ 

15.3 ± 1.5 

16.2 ±2.1 

13.7±0.7 

16.2 ± (1.6) 

(15.99 ± 1.18) 

(10.076 ±0.47) 

1 

+ 

A 

+ 

0.5 ±0.1 

0.4 ±0.1 

0.3 ±0.1 

0.6 ±0.1 

0.39 ±0.09 

0.204 ±0.040 

[K 7r+)g 

27.9 ±27.5 

34.2 ±25.7 

109.8 ±27.8 

— 

40.50 ±9.63 

— 

+ 

+ 

II 

bO 

1 

21.1 ±19.3 

5.1 ±5.6 

174.6 ±28.9 

— 

- 

15.5 ±2.8 

Background 

— 

— 

— 

17.2 ±5.3 

- 

— 

S-waves 

(81.7 ±0.8) 

(82.3 ±0.8) 

(81.2 ±0.7) 

(79.8 ± (12.0)) 

(83.23 ± 1.50) 

(97.1 ±3.9) 

Sum 

160.7 

134.5 

367.3 

66.3 

264.13 

122.9 


the {K and fit fractions individually. 

In Table V, we also list fit fractions from the previous analyses done by the experimental 
groups for a comparison with our result. Although each of the experimental groups ob¬ 
tained several models through their analyses, we do not try to exhaust all the models in 
the comparison. Rather we pick up some cases that are particularly interesting to compare 
with our results. Here, we tabulated three analyses results in Table V: the E791 Isobar 
model [5] on which our pseudodata are based, the FOCUS JC-matrix model analysis [8], 
and the CLEO QMIPWA [11]. We note that the definition for the fit fraction used by 
the E791 in Ref. [5] is different from what used here; the partial-wave amplitudes are not 
Bose-symmetrized in calculating the fit fractions, while our and the CLEO’s [11] fit fractions 
are from Bose-symmetrized amplitudes. However, we actually tabulated in Table V the fit 
fractions for the E791 Isobar model calculated by ourselves with our definition, and their 
errors are assumed to be the same as those given in Ref. [5]. In the original papers, the 
experimental groups presented the fit fractions of each of resonances considered. In order to 
(roughly) compare their results with ours, we sum the resonance contributions in the same 
partial wave incoherently, and the errors are summed in quadrature. The numbers enclosed 
by 0 are obtained by the incoherent (quadrature) sum. For the ' tt’*' fit fraction 

obtained by the incoherent sum, the A* (892) fit fraction dominates, and thus the interfer¬ 
ence effect would not be so large. A general trend in Table V is that all the analyses listed, 
and the E791 [5] and FOCUS [10] MIPWA results (not listed on the table) as well, are in 
fairly good agreement on the (A“7r+)p“^^^7r+, (A“7r+)^^'^^7r+, and “S-waves” fit fractions. 
Although the “S-waves” fit fraction from the CLEO QMIPWA is somewhat larger than the 
others, it is probably because the (7r+7r+);g"^A“ and Aq( 1430) fit fractions have been added 
incoherently. The other fit fractions are rather largely dependent on each of the analyses. 
For example, as we have mentioned, the FOCUS A-matrix model analysis gives rather large 
(A“7r’'')g~^'^^7r+ and (A“7r’'')g~^^^7r’'' fit fractions that interfere very destructively while the 
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TABLE VI. 

Bare ^ K 

decay fit fractions (%). 



Full 

Z 

Isobar 


45.0 ±9.2 

39.8 ±5.0 

— 

(A7r)p“^'^^7r+ 

13.1 ± 2.2 

5.8 ± 1.1 

— 

(iF7r)p=^/^7r+ 

0.3 ±0.2 

0.1 ±0.1 

— 

{TT~^7r^Yl^^^^ 

16.0 ±0.9 

42.7 ±0.7 

— 

(iF7r)J“^/^7r 

8.7 ± 11.7 

7.7 ±8.2 

— 

{ttttYs^'^K 

16.8 ± 19.1 

3.9 ±4.4 

— 


interference is moderate in the Full and Z models. 

With our coupled-channel analysis, it is interesting to examine bare fit fractions defined 
as follows. We first calculate a bare {ab)\^c partial width using the decay amplitude in which 
all contributions from the rescattering processes [the second term of Eq. (32); Fig. 6(b)] 
are omitted. In our coupled-channel description of the D'^ —)■ K~7r~^'ir~^ decay, we consider 
not only K~7r~^7r~^ but also K^7r^7r~^ in the hadronic intermediate states. Therefore, when 
calculating the bare partial width, we consider both of the states in the final state sum. 
(Precisely speaking, K^rj'7r~^, K^K^K~^, and the effective inelastic channels in the ttK p- 
and d-waves are also included in the intermediate states; their partial widths are rather 
small.) Then, the bare fit fraction for a given {ab)\^c is defined by the bare {ab)\^c partial 
width divided by the sum of the bare partial widths for all the considered {ahYi^c. Thus the 
sum of all the bare ht fractions is 100% by definition. The result is presented in Table VI. 
We leave the column for the Isobar model blank. This is because the —)■ TZc vertices 

of the Isobar model implicitly includes the rescattering effect, and we cannot eliminate the 
effect to extract the bare vertices. 

A remarkable feature in Table VI is the large fit fraction of in which /o(770) 

plays a major role. This is particularly true for the Z model. This fit fraction did not 
appear in Table V because {ti ^can appear only in intermediate states of the 
—)■ decay, and its contribution is genuinely a coupled-channel effect. The 

previous —)■ analyses cannot see this effect because they did not explicitly 

consider three-body dynamics. In order to see the contribution of the 

—)■ decay more clearly, we show in Fig. 14 the {7r~^7r~^) squared invariant 

mass spectrum of the Full model but couplings to are turned off. The figures 

clearly show the large contribution from {7r~^7rYp"^K^■ It increases the decay rate by ~7%, 
and significantly changes the shape of the spectra. We found a further enhanced contribution 
from {7r~^7rYp"^K^ for the Z model; it increases the decay rate by ~ 30%. Actually, the 
dominance of the {7r~^7r^)Y"^K^ ht fraction (~85%) was also found in a recent analysis of 
D~^ A'gTT+TT^ done by the BESIII Collaboration [12]. The JF“7r+7r+ and —)■ 

iF^TT+TT^ decays share the same hadronic dynamics, except for additional but much smaller 
doubly Cabibbo suppressed channels in the latter. Therefore, the large bare ht fraction of 
{7r~^7r^)Y"^K^ for the —)■ K~7r~^7r~^ decay seems natural in order to understand the decay 

mechanisms for —)■ K~7r~^7r~^ and D~^ —)■ Kg7r~^7r^ consistently. However, we have seen 
the rather large model dependence of the contribution from {7r~^7r^)Y"^. This would be 
largely due to the fact that we analyzed the D'^ —)■ data that do not contain this 

partial wave in the hnal state. Because the —)■ Kg7r~^7r^ decay contains (7r+7r°))T^.A° in 

the hnal state, a combined analysis of these two decay modes would signihcantly reduce the 
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FIG. 14. (Color online) Contributions of to the K~tt^ squared invariant 

mass spectrum in the left [right] panel. The Full model is shown by the red solid curve. The blue 
dashed curve is also from the Full model but couplings to are turned off. 




FIG. 15. (Color online) Contributions of the three-meson force (Fig. 4) to the K~7r^ [vr+Tr"’'] 
squared invariant mass spectrum in the left [right] panel. The Full model is shown by the red solid 
curve. The blue dashed curve is also from the Full model but the three-meson force is turned off. 
The magenta dot-dashed curve is from the Full model with all the rescattering effects turned off. 


uncertainty associated with (7r’''7r°)p 

Now let us study how much the three-meson force contributes to the K~7r~^7r^ 

decay. In Fig. 15, we compare the K~7t^ (tt+tt’'') squared invariant mass spectrum of the 
Full model with those from the same model but the three-meson force being turned off. The 
three-meson force suppresses the decay width by ~22%, and change the spectrum shape 
signihcantly, as seen in Fig. 15. The K~7i~^ spectrum is suppressed by the three-meson force 
at the JF*(892) peak region, and the 7r+7r+ spectrum is suppressed at higher We 

found that the effect of the diagram in Fig. 4(a) connected to {7r~^7rYp"^K^ is fhe most 
important among the three-body-type diagrams we consider. In the same hgure, we also 
show the spectrum from the Full model with all the rescattering processes [Fig. 6(b); second 
term in Eq. (32)] being turned off. Effects of the rescattering mechanisms are quite large; 
the decay width gets almost triplicated by the rescattering effect. The spectra are rather 
different from the blue dashed curves in Fig. 14 where {7r~^7i'^)Y"^K^ is turned off. Therefore, 
the hadronic rescattering through partial waves other than {7r~^7i'^)Y"^K^ also gives a major 
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FIG. 16. (Color online) Phase (upper) and modulus (lower) of —)• dressed vertex: 

(Left) TZf’^^c = i?5^7r+; (Middle) (Right) TZf’^^c = i?f^7r+. The red solid bands 

are from the Full model, the blue cross-hatched bands from the Z model, and the green bordered 
bands from the Isobar model; the band widths represents the errors. The unit of the modulus is 
arbitrary, but the relative magnitudes between the different bands in the three panels are given by 
the model. 


contribution. 

In a Dalitz plot analysis with an isobar-type model, it is assumed that a M* TZc [M*: 
heavy meson, meson resonance, etc.) decay vertex implicitly contains effects of rescattering 
mechanisms that are simulated by a complex coupling constant for the vertex. It is inter¬ 
esting to examine the extent to which this assumption is valid. In Fig. 16, the upper [lower] 
panels give the phase [modulus] of dressed —)■ TZc vertices defined in Eq. (32) for the Full 

model (red solid bands) and the Z model (blue cross-hatched bands), as a function of the 
momentum of the unpaired meson, pc- For comparison, we also show in the same hgure the 
D~^ TZc vertices from the Isobar model by the green bordered bands. The left, middle, 
and right panels are for TZf'^^c = TZf'^^c = and TZf'‘^^c = respectively. 

works as a seed to develop R and Kq( 1430) resonances, while and will develop 

.^*(892) and .^2(1430) resonances, respectively. Even though all of the models have been 
htted to the same pseudodata of the Dalitz plot distribution, they are rather different in 
the —)■ TZc vertices. In particular, the signihcant phase variations as a function of Pc 

in the Full and Z models are purely due to the three-body hadronic dynamics required to 
satisfy the three-body unitarity. The constant phase assumed in the isobar-type models is 
not justihed from this viewpoint. 


IV. SUMMARY AND FUTURE PROSPECTS 


In this work, we have performed a coupled-channel analysis of the pseudodata for the 
—)■ Dalitz plot. The pseudodata are generated from the isobar model of 
























the E791 Collaboration [5], and are reasonably realistic. As far as we know, this is the 
hrst coupled-channel analysis of a realistic Dalitz plot distribution for a D-meson decay 
into a three-pseudoscalar-nieson state. We have demonstrated that it is indeed possible to 
analyze this kind of high-quality data within a coupled-channel framework, and found lots of 
interesting results that are beyond what can be obtained with the conventional isobar-type 
model analyses. Let us summarize below what we have done and found in this work. 

In our build-up approach, we started with developing a two-pseudoscalar-meson inter¬ 
action model. With a suitable combination of contact interactions and bare resonance- 
excitation mechanisms, our two-pseudoscalar-meson interaction model successfully describes 
empirical tttt and tiK scattering amplitudes of low partial waves that are needed to ana¬ 
lyze the —)■ Dalitz plot. Poles associated with tttt and ttK resonances have 

been extracted, and most of them are in agreement with the PDG listings. Then using the 
two-pseudoscalar-meson interactions as building blocks, we developed a three-pseudoscalar- 
meson interaction model that describes the FSI of the —)■ K~7r~^7r~^ decay. The main 

driving force for the three-pseudoscalar-meson scattering process is the Z-diagrams and the 
dressed 77-propagator that appear as a necessary consequence of the three-body unitarity. 
These mechanisms do not contain any adjustable parameters once the two-pseudoscalar- 
meson interactions have been fixed using the tttt and nK scattering data. In addition, we 
considered mechanisms that are beyond simple iterations of the two-body interactions, and 
thus may be called a three-meson force. We introduced the three-meson force to an impor¬ 
tant channel in the FSI of the —)■ K~7r^7r^ decay, i.e., the vector-pseudoscalar channels. 

Guided by the hidden local symmetry model [21] that incorporates vector and pseudoscalar 
mesons in a chiral Lagrangian, we derived the vector-pseudoscalar interactions that work as 
a three-meson force. 

In our analysis of the pseudodata for the —)■ K~7r~^7r~^ Dalitz plot distribution, we 

basically used three models: Full, Z, and Isobar models. In the models, we took account 
of all of the channels (partial waves) that had been considered in the previous analyses 
of the same process. This includes (.A7r)^~^^^7r {{7r7r)g"^K) that is considered only in the 
FOCUS [8] (CLEO [11]) analysis. In addition, in the Full and Z models, we also considered 
where p(770) plays a dominant role. This partial wave can contribute to the 
—)■ A'“7r’''7r’'' decay only through the rescattering, and thus this contribution is a pure 
coupled-channel effect. A distinct feature of our models is that all of the two-pseudoscalar- 
meson resonances are included as poles of the unitary scattering amplitudes that ht well the 
empirical tttt and tiK amplitudes. With an appropriately selected set of htting parameters 
associated with bare D+ —)■ TZh {h = vr, K) vertices, the Watson theorem is satished by 
the models when the rescattering is turned off. The three-body unitarity is also taken care 
of within the Full and Z models. Unlike the previous analyses, we did not include a flat 
background amplitude. 

Our models ht the pseudodata with a reasonable precision. As far as the y^-value is 
concerned, the Full and Z models are signihcantly better than the Isobar model. Although 
this may be partly due to the larger number of adjustable parameters in the hts, it would 
also be because important mechanisms are considered in the better models. Indeed, we 
showed that the Isobar model can be signihcantly improved by just introducing the hadronic 
rescattering (x^/d.o.f. is reduced by ~29%), keeping the number of the htting parameters 
unchanged [Z(without p) model]. The inclusion of the pK channel further reduces y^/d.o.f. 
by 47%, and thus the importance of this channel seems rather clear. On the other hand, 
the three-meson force does not improve although it gives a signihcant ehect on the FSI. 
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Another important point is that we achieved the good hts using two-pseudoscalar-meson 
amplitudes hxed by the empirical tttt and ttK amplitudes. This is in sharp contrast with the 
previous analyses where the p-wave ttK amplitude is given by a sum of the Breit-Wigner 
functions that does not necessarily satisfy the Watson theorem. Also, the previous analyses 
commonly adjusted parameters associated with A’q( 1430) in their hts, and obtained widths 
signihcantly narrower than those in the PDG listings. In our analysis, on the other hand, 
we were able to obtain the reasonable hts with A’q( 1430) whose width is close to those from 
the PDG. 

We examined the partial-wave amplitudes and found that the Full and Z models have 
similar amplitudes, while the Isobar model has signihcantly diherent amplitudes, particularly 
for and . The and {Tr~^7r~^)g"‘^K~ amplitudes 

have rather large errors due to the fact that the htting parameters associated with these 
partial waves cannot be precisely determined by the data used in this work. We compared 
the amplitudes from our models with the {K~7r^)s'n'~^ amplitude from the 

E791 MIPWA. We found a good agreement between the Full and Z models, and the E791 
MIPWA for Mk-tt+ ^ 1-5 GeV. We stress that the hadronic rescattering plays an essential 
role to bring the phases of our models into agreement with those from the E791 MIPWA. 
On the other hand, the Isobar model, that maintains the Watson theorem, does not have a 
freedom to change the phases in the elastic region, ending up with a rather different solution. 
The partial-wave amplitudes were used to calculate the fit fractions that were compared with 
those from the previous analyses done by the experimental groups. We found a fairly good 
agreement between the analyses on the and £t fractions. 

For the and fit fractions, however, we found 

a rather large dependence on each of the analyses and also large errors. Even so, their 
coherent sum turned out to be consistent among our three models and also the 
fit fractions from the E791 and FOCUS MIPWA within greatly reduced errors. 

With the coupled-channel framework, we were able to study the bare fit fractions for which 
all the hadronic rescattering effects are absent. This quantity could be useful to study the 
intrinsic quark-gluon substructure of the D-meson. Remarkably, that does 

not show up in the usual fit fraction, gives a large bare fit fraction. This result may appear 
a bit surprising. However, considering that the D+ —)■ iF“7r+7r+ and —)■ decays 

share the same hadronic dynamics to a large extent, this finding is actually consistent with 
the recent BESIII analysis [12] that found a dominant fit fraction (~85%) of 
in the —)■ decay. 

We further studied coupled-channel effects. We found that the D”*" —)■ decay 

width gets triplicated by the rescattering effects in the Full model. We also found that the 
phases of the dressed —)■ vertices have rather large dependence on the unpaired 

momentum. The phase variation is a consequence of the explicit treatment of the three- 
body dynamics. In the conventional isobar-type model analyses of heavy or excited meson 
[M*) decay into three light mesons, it is assumed that the rescattering effects are small 
and/or the phases of the M* —)■ TZc vertices can be approximated by constants. Clearly, our 
analysis indicates that neither of these assumptions are supported from this more microscopic 
viewpoint, at least for the D~^ —)■ K~Tr~^7T~^ decay. With the above results as a basis, we can 
rather clearly conclude that explicit treatment of the hadronic FSI is essential for extracting 
partial wave amplitudes from Dalitz plots for D~^ iF“7r+7r+ and probably also many other 

processes. 

In future, we will apply our coupled-channel model to a combined analysis of the D~^ —)■ 
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and —)■ decays. The strength of the coupled-channel framework 

is to describe different processes in a unihed manner. Because different aspects of the 
hadronic dynamics appear in different processes, the combined analysis is a powerful method 
to understand the hadron dynamics with smaller model-dependence. This is why a com¬ 
bined analysis has become standard in the study of the baryon spectroscopy. For example, 
ttA^, 'yN —)> ttN, rjN, KA, KH reaction data are analyzed with a coupled-channel model in 
a unihed manner to extract the nucleon resonance properties in Ref. [34]. This direction 
should also be pursued to better understand the hadronic dynamics in heavy-meson decays 
and meson resonances. Getting back to the future combined analysis of —> K~7r~^7r~^ and 

we expect to better understand the role of (7r+7r°)|r^.R° because contributions of 
this partial wave can be directly seen in the D~^ —)■ Dalitz plot. Although we found 

the very important contribution of {7r~^7r^)^jr^K^ to the D'^ —)■ K~7r~^'n'^ decay through the 
FSI, this hnding is based on the indirect information. Also in the combined analysis, we 
would be able to better extract the {K7r)g~^^‘^7r and {7r7r)g"‘^K amplitudes that have been de¬ 
termined with large uncertainties in this work. By simultaneously analyzing the D~^ decays 
with differently charged hnal states, it will be possible to better separate contributions from 
these different isospin states. It will be interesting to examine how the partial-wave ampli¬ 
tudes and ht fractions obtained in this work will be in the combined analysis. Finally, we 
will also study if a three-meson force plays an important role in understanding the hadronic 
FSI of the two Zi)+-meson decays in a unihed manner. 
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Appendix A: Three-meson force based on Hidden Local Symmetry Model 

In this Appendix, we hrst present a set of Lagrangians from the hidden local symmetry 
(HLS) model [21]. Then we present expressions for potentials, derived from the Lagrangians, 
between vector-mesons and pseudoscalar-mesons. These potentials work as a three-body- 
force in a three-pseudoscalar-meson system. Finally, we present the potentials in a partial- 
wave form that is useful for numerical calculations. 


1. Lagrangians from the HLS model 

We use symbols P and V to denote octet pseudoscalar-mesons and nonet vector-mesons, 
respectively. The pseudoscalar meson helds in the SU(3) matrix form are 
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where Aq is a Gell-Mann matrix, while the vector meson nonet is given by 


V 


1 8 
a=0 


1 

71 



p 

K*- 


V2^_ ^ 
K*o 


V2 


u 



(A2) 


where Aq = y |l (1: unit matrix), and the ideal mixing between the nentral vector mesons 

is assumed. When P and V are enclosed in the trace symbol, they are fields of the SU(3) 
matrix form. Otherwise, e.g., they are in brackets, they are understood to be one of particles 
contained in the SU(3) matrix elements. It is convenient to work with isospin states rather 
than the charge states used in Eqs. (Al) and (A2). For the relation between the two sets 
of the basis, we employ a convention where the charge states are the same as their isospin 
states (|//^)) with some exceptions that need additional phases as follows: 

|p+) = -|/ = 1, r = 1), \K-) = -\I = 1/2, 1‘ = -1/2), 

|7r+) =-|/= l,r = 1), =-|/= 1/2,/-'=-1/2) . (A3) 

In what follows, \P) and \V) are understood to be isospin states rather than charge states. 
Also, we use curly symbols to denote creation or annihilation operators. For example, V is 
the annihilation operator contained in the field P, and its normalization is (0|P|Pa) = Aa/2. 

The mesonic interaction Lagrangians we use are those from the HLS model [21, 35]. The 
VPP interactions are (with the Bjorken-Drell convention for the metric) 

Cvpp = 2igTT[Vf,{df^PP - Pd^P)] , (A4) 

where the trace is taken in the SU(3) space. The coupling g is related to the pmi coupling by 
g = gpiT-K, and we use gp^n = 6.0 determined from the p —?■ tttt decay width. The Yang-Mills 
type Lagrangian, from which we use VVV interactions, is 

^VM = -^Tr [Fp^F^ , (A5) 

with 

Fp. = dpV, - d,Vp - tg[Vp, W] . (A6) 

The VVP interactions are given by 

CvvP = p'Ge^^^'^Tr [d^V^d^VsP] , (A7) 

where we use the convention, = -|-1. The coupling is G = —3/(47r^/7r) with being 
the pion decay constant. In our numerical calculations, we use \g'^C/2\ = g^^T^p ~ 0.012 from 
an analysis by Durso [36] on the decay u ^ Tip ^ vry. 


2. YP —)• V'P' potentials 

We consider a process V{pv, ey, Jy, Iy)P{pp, Ip, Ip) —)■ V'{pv', ey', Iv>,Py,)P'{pp^,Ip',Pp,) 
where the variables in the parentheses specify each particle’s state such as the four- 
momentum (p), polarization (e), isospin (J) and its ^(-component {P)- The 0th component 
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of the four-momentum is related to the spatial part by = \/Px + where rrix is the mass 
for a particle x. The potential diagrammatically represented in Fig. 4(a) is derived from the 
Lagrangians in Eqs. (A4) and (AS) following the unitary transformation method [37], and 
is given by 


Id 


Fig.4(a) 




— m?/ 
\-i 


.(_l)A+2V+^w.-2/h+A'+i^(2/p/ l){2Iv + 1) 




I 


X {{pv + Pv) ■ {pp + PP')ev ■ A/' - (pv + q) ■ (*v'i.PP + Pp') ' ev' 
+ {q- Pv) ■ ey(pp + pp') ■ e*y>} , 


(AS) 


where q = Pv ~ Pv = Pp' ~ Pp^ and and ly^^ are the exchanged vector-meson mass 
and isospin, respectively. We have used the Racah and Clebsch-Gordan coefficients denoted 
by lR(jij 2 ' 7 j 3 ; J 12 J 23 ) and (jimij2P2-2|'7M), respectively. The propagator for the exchanged 
vector-meson is more explicitly written as 


2[(pv-py/)2 


m 


Ve= 


{Pp' -Pp)^ -m 


(A9) 


as specihed by the unitary transformation [37]. The effective coupling gvvvex is given by 

g Tr[(V)[(V'),(Ve.)]] 


9VVVe:, — 


\/2 {IvIv'^Ve:,IvJIvIv) 


(AlO) 


with (V) = (0|V|R), (V') = (R'|V|0) and (Vex) = (Rea;|V|0). The effective couplings gwv^^ 
as well as gv^xPP' and gvVe^P' appearing below are independent of isospin ^-components. 
9v,.pp' is given by 


9v^^pp' — 


g TV|(H,)|(P).(P)|] 

\/2 (I Vex ^Vex ^P^pl^P'^P') 


(All) 


with (Vex) = (0|V|Wx), {V) = {0\r\P) and {V) = {P'\V\0). 

Another potential diagrammatically represented by Fig. 4(b) is derived from the La- 
grangian in Eq. (A7) and is given by 




Fig.4(b) 


9VVexP'9VVexP 
9^ - ^lex 


\/ ( 2 /y + l)( 2 /y/ -|- 1 ) 


X 


W{IyJpIpJy-Iy^J){IyFyIpFp\IF){Iy,Py,Ip,Fp,\IP 


X {pyi ■ pp'{e*y, ■ Pvev ■ Pp - e*y' ■ evPv ■ Pp) + Pv ■ Pvi^y, ■ evpp ■ pp' - ep, ■ pp^ev ■ pp) 
+ey • Pv'(e*v' ■ Pp’Pv • Pp - ' PvPp ■ Pp')} , (A12) 

where the propagator for the exchanged vector-meson is 

1 1 r 1 1 ^ 


4 

y ez 


+ 


{pv - Pp )2 - 


9^ - ^Vex 2 [{pv- Pp')^ - m 
The effective coupling gvVexP' is given by 

'IV1(P')({V>(UJ + (U,){V»1 


gvv^^p' — 2c/< 


LdTTp 




(A13) 


(AM) 
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with {V) = (P'|P|0), (V) = (0|V|i/) and {Vex) = (Ka;|V|0). gvVe^P is obtained by exchang¬ 
ing labels V ^ V and P' —)■ P in the rhs of Eq. (A14), and the meaning of the matrix 
elements are (V') = (E'|V|0), (P) = (0|P|P) and {Vex) = (0|V|14x). 

Following Ref. [25], we mnltiply the following form factor to the potentials of Eqs. (AS) 
and (A12): 

(A15) 

where Asmf is a cutoff to be determined by fitting data. We use the same cutoff value for 
all the potentials of Eqs. (AS) and (A12). We checked numerical values of the potentials by 
comparing our calculation with Fig. S of Ref. [25]. 


3. Partial-wave expansion 

In order to implement the VP —)■ VP' potentials of Eqs. (AS) and (A12) into the 
scattering equation of Eq. (17), we need to expand the potentials in terms of the partial- 
wave representation as, 

{V{py>, ev, Iv, Iv')P\Vp'. Ip', Ip')\VFi,A\V{pv, ev, ly, Iv)P{Pp, Ip, Ip)) 

= Y. {IvIvIpIp\TPy + Pp){Iv>Py,IpP^p,\TPy + Pp) 

TJJHPl'l'^ 

where JPT are the total angular momentum, parity, and total isospin, respectively, and I 
[I') is the orbital angular momentum of the relative motion of RP {VP'). Inverting this 
equation, we obtain 


VV'i7pv)fiPv'l\Pv\) = (IvHJpmTT’)(Ivd^,Ip,I’i,.\Tr) 

JZ TZ JZ JZ 

X (-1)^+^ +^v+ew (le^ J - ev|/0) (ley/J - ey\l'eyi - ey) J 

X {V{py>, ey,, ly,, Py,)P'{pp,, Ip,, Pp,)\Vf,,a\V{ py , ey, ly, I^)P{pp, Ip, Pp)) , (A17) 

where py is taken along the z-axis. This partial-wave form of the potentials is plugged in 
Eq. (26). 


47r 


2 /' + 1 


Appendix B: Model parameters 
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TABLE VII. Model parameters for the -kK partial-wave scattering with the angular momentum L 
and the isospin I. The parameters are defined in Eqs. (1), (2), ( 8 ), and (10). Eor each partial wave 
specified by {L,/}, masses couplings {ghih 2 ,Ri)^ cutoffs {ch^h 2 ,Ri) bare R 

states, Ri, and hi and /12 are particles in a two-pseudoscalar-meson channel. Couplings {hh^h 2 Mh 2 ) 
and cutoffs ( 6 / 11 / 12 ) are for the contact interactions. Masses (mi, m 2 ) are for two “particles” in 
the effective inelastic channel, except for the {L,/}={ 0 ,l/ 2 } wave for which (mi, m 2 )={mK, rrirf'). 
The superscripts, LI, of the parameters are suppressed for simplicity. The masses and cutoffs are 
given in the unit of MeV. The hyphens indicate unused parameters. 


R{L,I} 

K*o {0, 1/2} 

{0, 3/2} 

Kt {1, 1/2} 

K*2 {2, 1/2} 


1391 

— 

1081 

3070 

9ttK,Ri 

-6.28 

— 

0.52 

0.18 

^■kK,Ri 

609 

— 

1973 

1954 

9h\h2iR\ 

4.30 

— 

-0.00 

0.97 

Chih2,Rl 

1966 

— 

1706 

1035 

rriR^ 

1767 

— 

1580 

— 

9-kK,R2 

8.22 

— 

1.84 

— 

^ttK,R2 

395 

— 

395 

— 

9hih2,R2 

-4.87 

— 

-3.00 

— 

^hih2,R2 

395 

— 

411 

— 

rriR^ 

— 

— 

1750 

— 

97tK,R3 

— 

— 

0.23 

— 

^tvK,R3 

— 

— 

1316 

— 

9hih2,R3 

— 

— 

2.65 

— 

Chi/i2,/?3 

— 

— 

395 

— 

^ttK ,7vK 

1.31 

0.45 

— 

— 

Kk 

710 

1973 

— 

— 

mi 

494 

— 

591 

100 

m 2 

958 

— 

662 

1049 


35 







TABLE VIII. Model parameters for vrvr partial-wave scatterings. The coupling h' is defined in 
Eq. (3). The other features are the same as those in Table VII. 


R{L,I} 

/o {0, 0} 

{0, 2} 

P{1, 1} 

/2 {2, 0} 

rriR^ 

1166 

— 

850 

1561 

9'K'k^R\ 

5.97 

— 

1.02 

-0.32 

Ctttt,/?! 

1162 

— 

805 

962 

9kk,Ri 

-2.19 

— 

-0.18 

0.19 

^KK,Ri 

1973 

— 

395 

1216 

rriR^ 

1627 

— 

1551 

— 

9'inT,R2 

-5.23 

— 

0.49 

— 

^7nT,R2 

1973 

— 

1973 

— 

9kk,R2 

11.99 

— 

3.74 

— 

^KK,R2 

533 

— 

395 

— 

^TTTTjTTTr 

0.47 

0.11 

— 

— 

h' 

— 

0.21 

— 

— 

^TTTT 

897 

913 

— 

— 


36 







TABLE IX. Parameters determined by fitting D'^ —^ K Dalitz plot pseudodata; —)• 

KTZ^ ’ {h = TT, K) bare coupling i/, 2 /) and phase ( 0 ^^z., 2 /) in degrees, as defined in Eq. ^33). 

stands for Lth bare R state with the spin L and isospin I. stands for the spurious state 
(see Sec. IIB for the definition) that decays into two pseudoscalar mesons {ab) with the orbital 
angular momentum L and total isospin I. The total spin, parity, total isospin, and hTZ relative 
orbital angular momentum are J=0, P=+l, T=3j2, and l=L, respectively for all the parameters, 
and thus labels of JPT and I are suppressed in the table. The second, third, fourth, and fifth 
columns show the parameters for the Eull, Z, Isobar, and Z(without p) models, respectively. The 
hyphens indicate unused parameters. The cutoff for the three-meson force in the unit of MeV is 
denoted by Asmf [see Eq. (A15) for definition]. 



Full 

Z 

Isobar 

Z (without p) 

'/’ttKOI 

-72. ±5. 

-80. ±4. 

2. ±3. 

-24. ± 6. 


10.54 ±0.57 

24.32 ± 2.94 

10.83 ±0.60 

7.81 ±0.79 

CttROI 

0.00 ± 0.28 

0.00 ± 0.87 

8.37 ±0.59 

12.34 ± 1.16 


0.00 ± 0.01 

0.00 ± 0.03 

0.00 ± 0.00 

0.00 ± 0.00 

‘(’Trfill 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 


1 (fixed) 

1 (fixed) 

1 (fixed) 

1 (fixed) 


1.48 ±0.31 

4.08 ± 0.83 

0.03 ±0.26 

0.19 ±0.24 


1.34 ±0.46 

6.55 ± 1.28 

0.03 ± 0.43 

0.39 ±0.50 

4>nR'2^ 

82. ± 13. 

136. ± 24. 

-131. ±9. 

129. ±6. 

C'TTijf 

0.27 ±0.08 

0.37 ±0.10 

0.21 ±0.04 

0.34 ±0.04 

4>kr^^ 

-44. ± 8. 

178. ± 7. 

— 

— 

Ckr\^ 

0.00 ± 0.00 

19.37 ±3.12 

— 

— 

CrrI^ 

17.27 ± 1.48 

0.00 ± 0.84 

— 

— 


121. ± 19. 

62. ± 16. 

-133. ±4. 

-170. ±8. 


0.90 ±0.52 

2.06 ± 1.12 

2.61 ±0.26 

3.00 ±0.36 


173. ± 30. 

36. ± 42. 

-78. ± 5. 

-114. ±7. 


1.78 ± 0.79 

2.09 ± 1.20 

4.67 ±0.30 

4.97 ±0.48 

Asmf 

2563. ± 722. 

— 

— 

— 
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